PROCEEDINGS 

°F SCIENCE 



O 
O 
(N 

> 
O 



> 

oo 



o 

o 
U 



X 



"Oscillating" components in the BL Lac object 
0716+714? 



S. Britzen™, V. Meyer'', A. WitzeP, I. Agudo^ M.F. Aller^ H.D. Aller' , A. Eckarr 
J.A. Zensus'' 

"Max-Planck-InstitutfurRadioastronomie,AufdemHilgel69, 53121 Bonn, Germany 

^Astronomy Department, University of Michigan, Ann Arbor, MI 48109-1042, USA 

'^I. Physikalisches Institut Koln, Universitdt zu Koln, Ziilpicher Str. 77, 50937 Koln, Germany 



and 



awit zel@mpif r-bonn . mpg . de iagudoSmpif r-bonn . mpg . de mf a@umich . edu 



E-mail: sbr it zen@mpif r-bonn . mpg . de vmeyergmpif r-bonn . mpg . de, 



haller@umich . edu eckartSphl . uni-koeln . de 
azensus@mpf ir-bonn . mpg . de 



15 VLBA observations of the BL Lac object 0716+714, performed between 1994.67 and 2006.40 
at 15 GHz, have been analyzed. Part of the data (12 epochs) were obtained within the Mojave- 
survey project and reanalyzed by us. We present a new motion scenario for jet model component 
motion in this source which suggests that no longterm-outward motion but rather an oscillation of 
components around an average core separation is taking place. Although no significant outward 
motion in the core separation can be found, motion with regard to the position angle is observed. 
We give lower limits for the derived apparent motions and Doppler factors and compare them with 
values published in the literature. We find a relation between the total flux-density evolution at 
14.5 GHz (UMRAO data) and the position angle changes observed for the jet component closest to 
the core. We suggest a significant geometric contribution to the longterm flux-density variability 
in 0716+714. 
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1. Introduction 

The BL Lac Object 0716+714 (redshift > 0.3, [pO|]) shows flux-density variability at all wave- 
length ranges that have been investigated so far (e.g., [[T9|]). It is still the only AGN for which 
a simultaneous change in the mode of the radio/optical variability has been observed on intiaday 
time scales ([p^]) and is thus so far the best candidate for an intrinsic origin of intraday variability 
(IDV, [ p^ and references therein). 0716-1-714 has been investigated with VLBI at various radio 
frequencies with different arrays. On parsec-scales it reveals a North-South directed jet with sev- 
eral components imbedded (e.g., [^). On kpc-scales, the jet is more East-West oriented with some 
evidence for halo-emission (e.g., ^). A rather large range in proper motions (0.05 mas/yr - 1.1 
mas/yr) has been reported by different authors based on the investigation of an increasing core 
separation with time modelled for jet components in this source (e.g., [||], [^, [[l8|], [|l0|], [pl|], 
[1T6|], [p^, [Q]). Based on these results it is not clear, whether 0716-1-714 reveals the slow velocities 



typical for BL Lac objects (e.g., []g]) or faster velocities more typical for quasars (e.g., []12|]). We 
here present first results of an (re-)analysis of 15 epochs of VLB A observations revealing evidence 
for a new motion scenario in the jet of this source. In this paper, we concentrate on the description 
of the results based on the 15 GHz data. A paper describing these results in more detail and the 
implications of a multifrequency analysis is in preparation. In contrast to previous studies - that 
relied mainly on investigating the core separation evolution - we study the position angle evolution 
with time as well. Again in contrast to previous studies, based on our new identification scenario, 
we find no evidence for any longterm-outward motion but instead find that all modelled jet com- 
ponents remain at similar core separations and move with regard to the position angle. We call this 
kind of motion "oscillatory". A similar motion phenomenon has already been observed in the jet 
of 1803+784 ([||,[|l|],[01). 

2. The observations 

0716+714 has been observed in 15 VLBA observations at 15 GHz between 1994.67 and 
2006.40. 12 of these images have been obtained by the MOJAVE / 2 cm- VLBA survey group 



(MOJAVE, e.g., [14], 2cm Survey, e.g., [12]). 0716+714 was used as calibrator source for ob- 
servations targeted on NRAO 150 ([|l|]) in the epochs 2004.64, 2004.97, and 2005.25. Part of the 
data presented here have already been analyzed by [Q]. We show the results of a reanalysis of all 
the data sets. Modelfits within the difmap package (e.g., [[I5|]) have been performed starting from 
a point-like model in all the observations. The final number of jet components necessary to fit 
the data adequately was reached when adding another jet component did not lead to a significant 
improvement of the value of chi-square. The uncertainties of the model component parameters 
have been determined by comparing the parameter ranges obtained by performing modelfits with 
different numbers of model components (±1 component). 

2.1 Component identification 

The component identification was done based on the assumptions that changes of the flux- 
density, the core separation, and the position angle of the modelled jet components should be small 
on time scales between adjacent epochs. We show an example for the component identification in 
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Fig.|l] (left). This identification across the epochs leads to a new motion scenario (described in the 
next section) for this source compared with earlier scenarios published in the literature. 



Clean LL map. Array: BFHKLMN0P5 
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Figure 1: We show the source structure in a map obtained at 15 GHz with the contours of the model 
superimposed (left). The inner three jet components are labelled. The core separation as function of time 
based on the 15 GHz data is shown in the figure on the right. A linear regression based on a previous 
identification scenario (for a subsample of the data) is shown with a dashed line (components c2-cl2, [Q|). 
Components a-x result from the new modelfits (red). The linear regression for this new identification is 
plotted as solid line. The absolute values of the apparent motion with uncertainties in brackets are listed. 



3. Results 

3.1 Oscillatory components 

In Fig.|l] (right) we show the core separation as a function of time for all model-fitted com- 
ponents from the 15 GHz data. In addition to our own results, the results published by [Q] for a 
subset of the data presented here are shown as well. The component identification across the ob- 
serving epochs by Bach et al. is based on a combined data set at 5, 8, 15 and 22 GHz. For a better 
comparison with our component identification, we here only show the 15 GHz data with the linear 
regression made by [Q] (dashed line) superimposed. The larger amount of data analyzed here leads 
to a different component identification (solid line). The fast apparent velocities found by [Q] (noted 
in Fig.|I| (right)) can not be confirmed within our new identification of model components. Instead, 
we find that our derived values for the apparent proper motions of model components are consistent 
with no longterm outward-motion if we calculate these values based on the core separation versus 
time. All jet components show significant changes with regard to the position angle. We show this 
relation for the inner three jet components in Fig]^ (left). In Fig.^ (right) we show the motion of 
these model components in rectangular coordinates. 
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Figure 2: The position angle as function of time is shown for the inner three model components (a - c) (left). 
On the right, the inner component's motion are shown in rectangular coordinates (5-point running average 
for a, 4-point running averages for b and c). 

3.2 Apparent velocities and Doppler factors 

Some of the model components reveal apparent inward motion, which is not uncommon in 
AGN. However, none of these values is of 3 sigma significance. For the apparent velocities derived 
for the motion based on the core separation evolution with time only, we find values between 
O.libO.lc and 1.4ib2.2c. The Doppler factors for this motion are between 1.1 it 1.2 and 3.0±4.4 
depending on the assumed angle to the line of sight of 4.9° or ~0°, respectively (angles to the 
line of sight taken from [Q]). For the non-radial motion, we derive preliminary lower limits for the 
apparent velocities between 5±3c and 10±5c and Doppler factors between 8ib4 (iimer component 
a) and 11±5 (outer component d) (4.9°) and 9±5 (a) and 20±10 (d) (~ 0°) for the individual 
components. These Doppler factors have been derived taken the motion in rectangular coordinates 
into account. They compare well with Doppler factors derived based on e.g., Inter-Day variability 

3.3 Correlation between total flux-density evolution & position angle evolution 

In Fig.^ we show a composition of several data sets: the total flux density at 14.5 GHz (black 
dots, UMRAO data), the flux-density of the core (dot-dashed), the flux-density of component a 
(closest component to the core, dashed), and the position angle for component a (dotted) as function 
of time.We find that most of the total flux-density is contained in the core component. We also 
find evidence for a positive relation between the total flux-density evolution with time and the 
position angle evolution of component a with time. This might suggest that the longterm flux- 
density evolution in 0716-1-714 may have a non-negligible geometric contribution. 

4. Conclusions 

We present a new model component motion scenario for 0716-1-714 with the following impli- 
cations: 

• In this scenario the jet components show no 3 sigma evidence for longterm-outward motion. 
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Figure 3: The total flux-density at 14.5 GHz (UMRAO data) is shown as a function of time. The VLBI 
core flux-density (15 GHz, dot-dashed line) and the VLBI flux-density of component a (dashed Une) are 
superimposed. The left y-axis measures the flux-density. In addition, the position angle for component a 
(component closest to the core) is shown as a function of time (dotted line). The right y-axis denotes for the 
position angle. 



instead they remain at similar core separations and seem to change significantly in position angle 
within our computed uncertainties. 

• Preliminary lower limits for the apparent velocities (with regard to the position angle motion) lie 
between 5c (4.9°) and 10c (~0°), the Doppler-factors between 8-11 (4.9°) and 9-20 (~0°) for the 
iimer to outer components. These values compare well with the limits derived e.g., from variabiUty 
arguments. 

• We find a relation between the total flux-density and the position angle of the component closest 
to the core (a), which might suggest a non negligible geometric origin of the longterm flux-density 
variability in 0716-1-714. 

A model to explain the origin of this correlation (based on rotation or precession) is currently in 
preparation. In addition, we are investigating a large set of multifrequency data (5, 8, 22, 43 GHz) 
to study the frequency dependence of the presented relations. 

The usually observed outward motion of jet components in AGN is - according to our analysis - 
not observable in the BL Lac object 0716-1-714. Since BL Lac objects most Ukely constitute the 
class of AGN with a jet directed almost in the direction of the line of sight, projection effects, e.g. 
the amplification of curved structures due to the small angle to the Une of sight, might play an 
important role and will have to be disentangled before we draw final conclusions. 
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